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Abstract 

Resistive  companion  (RC)  modeling  is  an  easy-to-use  approach  for  electric  circuit  simulations.  With  a  RC  numerical  solver,  simulations 
of  complex  electric  systems  can  be  achieved  based  on  RC  models.  In  this  study,  the  construction  of  RC  battery  models  is  investigated.  A 
general  battery  model  and  a  nickel-metal  hydride  cell  model  have  been  built.  Simulations  of  RC  battery  models  on  virtual  test  bed  (VTB) 
are  presented  and  analyzed.  It  is  shown  that  RC  modeling  provides  a  flexible  and  powerful  way  for  the  simulation  of  battery  systems. 
©  2001  Elsevier  Science  B.V.  All  rights  reserved. 
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1.  Introduction 

Battery  becomes  more  and  more  important  in  our  life  with 
the  growing  popularity  of  portable  electric  devices,  e.g. 
camera,  cellular  phone,  and  laptop,  etc.  Some  high-energy 
batteries  (e.g.  lithium  ion  battery  and  nickel-metal  hydride 
battery)  are  under  active  research  at  present,  which  hope¬ 
fully  will  lead  to  some  breakthroughs  in  the  power  storage 
techniques  in  the  future.  Due  to  the  complex  chemical  and 
physical  processes  involved,  the  behavior  of  battery  is 
usually  harder  to  predict  compared  with  that  of  electric 
and  mechanic  devices  [1].  A  widely  used  battery  model  is 
described  by  the  following  equation 

E  —  Eo  +  IR[nt  (1) 

In  the  above  equation,  it  is  assumed  that  a  battery  has  a 
constant  open-circuit  potential  and  the  potential  of  a  loaded 
battery  varies  linearly  with  the  applied  current.  This  model  is 
actually  oversimplified,  i.e.  it  cannot  reflect  many  charac¬ 
teristics  of  a  real  battery,  such  as  only  limited  amount  of 
energy  is  available  from  a  battery,  and  the  open-circuit 
potential  of  a  battery  changes  with  the  state  of  charge 
(SOC),  etc.  Hence,  it  cannot  be  used  in  serious  simulations. 

First-principle  battery  modeling  becomes  popular  in 
recent  years.  It  is  an  effective  approach  for  battery  design 
and  optimization.  Since  detailed  chemical  and  physical 
processes  are  considered,  first-principle  battery  models 
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are  generally  complicated  (with  a  large  set  of  differential/ 
algebraic  equations  and  tens  or  hundreds  of  parameters).  As 
a  result,  these  models  require  significant  computational 
power  and  cannot  be  integrated  easily  with  other  simulation 
programs  to  investigate  the  behavior  of  whole  power  sys¬ 
tems.  In  fact,  most  first-principle  battery  models  in  literature 
are  coded  in  Fortran  for  standalone  simulations  and  are 
usable  for  only  limited  operations. 

Circuit  simulation  has  been  widely  used  in  the  electrical 
engineering  held.  There  are  many  circuit  simulation  pro¬ 
grams  available  today,  e.g.  widely  used  SPICE  [2].  Recently, 
a  heavily  funded  project,  virtual  test  bed  (VTB),  provides 
another  choice  of  electric  circuit  simulation  programs 
(downloadable  from  http://web.engr.sc.edu/vtb/).  VTB 
was  originally  designed  to  simulate  an  electric  ship.  Later 
on,  it  was  expanded  to  provide  more  capabilities  to  simulate 
general  engineering  systems,  e.g.  hybrid  electric  vehicles. 
VTB  has  two  major  components:  design  environment  and 
visualization  engine.  The  first  component  (see  Fig.  1)  pro¬ 
vides  a  visual  design  platform  and  performs  numerical 
calculations.  The  second  component  (see  Fig.  2)  accepts 
the  results  from  the  first  component  and  presents  them  in 
dynamic  plots.  A  comprehensive  library  of  models  has  been 
developed  for  VTB,  which  permits  simulation  of  many 
engineering  systems.  Detailed  information  about  the  VTB 
project  can  be  found  on  the  VTB  web  site. 

After  a  system  has  been  designed  in  the  VTB  environ¬ 
ment,  the  numerical  solver  of  VTB  utilizes  the  information 
contained  in  the  worksheet  and  the  resistive  companion  (RC) 
models  of  the  relevant  devices  to  predict  the  dynamic 
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Nomenclature 

Clneg 

specific  surface  area  of  the  negative  elec¬ 
trode  (cm2/cm3) 

Opos 

specific  surface  area  of  the  positive  elec¬ 
trode  (cm2/cm3) 

A 

^neg 

geometry  area  of  the  negative  electrode  (cm2) 

A 

^pos 

geometry  area  of  the  positive  electrode  (cm2) 

b 

equivalent  current  source  of  an  electric 
device  (A) 

B 

current  source  vector  (A) 

concentration  of  KOH  electrolyte  (mol/ 
cm3) 

^-e,ref 

reference  concentration  of  KOH  electrolyte 
(mol/cm3) 

CMH 

concentration  of  hydrogen  in  metal  hydride 
(mol/cm3) 

^MH,max 

maximum  concentration  of  hydrogen  in 
metal  hydride  (mol/cm3) 

cMH,ref 

reference  concentration  of  hydrogen  in 
metal  hydride  (mol/cm3) 

CNi(OH)2 

concentration  of  Ni(OH)2  (mol/cm3) 

^Ni(OH)2,max 

maximum  concentration  of  Ni(OH)2  (mol/ 
cm3) 

cNi(OH)2,ref 

reference  concentration  of  Ni(OH)2  (mol/ 
cm3) 

Chat 

capacity  of  the  battery  (A  h) 

E 

potential  of  a  loaded  battery  (V) 

Eo 

standard  potential  of  a  battery  (V) 

Ee  q 

open-circuit  potential  of  a  battery  at  a 
particular  state  of  charge  (V) 

/ 

linear  or  nonlinear  equations 

F 

Faradic  constant  (98467  C/eq) 

F 

governing  equations  of  voltage  variables 

8 

equivalent  conductance  of  an  electric  de¬ 
vice  (S) 

G 

conductance  matrix  (S) 

i 

current  through  an  electric  device  (A) 

i 

vector  of  current  (A) 

4, ref 

rate  of  oxygen  reaction  at  the  reference 
state  (A/cm2) 

Cell 

current  through  a  cell  (A) 

4 

current  on  terminal  k  (A) 

io,k 

exchange  current  density  of  reaction  k  (A/ 
cm2) 

ji 

current  density  due  to  the  nickel  reaction  on 
the  positive  electrode  (A/cm2) 

h 

current  density  due  to  the  oxygen  reaction 
on  the  positive  electrode  (A/cm2) 

h 

current  density  due  to  the  hydrogen  reaction 
on  the  negative  electrode  (A/cm2) 

h 

current  density  due  to  the  oxygen  reaction 
on  the  negative  electrode  (A/cm2) 

J 

Jacobian  matrix  with  respect  to  voltage 
variables 

4cg 

thickness  of  the  metal  hydride  electrode 
(cm) 

^pos 

thickness  of  the  nickel  electrode  (cm) 

ly,  neg 

equivalent  thickness  of  metal  hydride 
material  (cm) 

^y,pos 

equivalent  thickness  of  nickel  active  mate¬ 
rial  (cm) 

L 

lower  triangular  matrix 

Emu 

loading  of  metal  hydride  material  (g/cm2) 

LNi(OH)2 

loading  of  nickel  active  material  (g/cm2) 

m 

temporary  vector  in  matrix  manipulation 

44/ni 

molecular  weight  of  metal  hydride  material 
(g/mol) 

4fNi(OH)2 

molecular  weight  of  nickel  active  material 
(g/mol) 

n 

electron  transfer  number  of  the  whole 
battery  reaction 

Po2 

pressure  of  oxygen  (atm) 

P02  ,ref 

reference  pressure  of  oxygen  (atm) 

R 

ideal  gas  constant  (8.3143  J/mol/K) 

E'mt 

internal  resistance  of  a  battery  (Q) 

soc 

state  of  charge  of  the  battery 

t 

independent  time  variable  (s) 

T 

temperature  (K) 

u 

vector  of  independent  control  variables 

U 

upper  triangular  matrix 

64q,A:,ref 

equilibrium  potential  of  reaction  k  at 
reference  reactant  concentrations  (V) 

V 

voltage  across  an  electric  device  (V) 

V 

vector  of  voltages  (V) 

^cell 

voltage  through  a  cell  (V) 

Vk 

voltage  on  terminal  k  (V) 

Fgas 

gas  volume  in  a  nickel-metal  hydride  cell 
(cm3) 

y 

vector  of  internal  state  variables 

Greek  letters 

A/cell 

change  of  current  due  to  the  perturbation  of 
the  voltage  (A) 

A’Tcell 

perturbation  of  the  voltage  (V) 

^^neg 

potential  different  at  the  solid-liquid  inter¬ 
face  on  the  negative  electrode  (V) 

^^pos 

potential  different  at  the  solid-liquid  inter¬ 
face  on  the  positive  electrode  (V) 

^neg,l 

potential  in  the  liquid  phase  on  the  negative 
electrode  (V) 

^neg,s 

potential  in  the  solid  phase  on  the  negative 
electrode  (V) 

^pos,l 

potential  in  the  liquid  phase  on  the  positive 
electrode  (V) 

^pos,s 

potential  in  the  solid  phase  on  the  positive 
electrode  (V) 

P  MH 

density  of  metal  hydride  active  material  (g/ 
cm3) 

PNi(OH)2 

density  of  nickel  active  material  (g/cm3) 
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Fig.  1.  The  VTB  design  environment. 
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Fig.  2.  The  VTB  visualization  engine. 
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Fig.  3.  RC  black  box  model  of  an  electric  device. 


behavior  of  the  system.  For  each  electric  device,  a  black  box 
RC  model  with  some  external  terminals  (see  Fig.  3)  [3]  can 
be  used  to  describe  its  electric  behavior.  The  governing 
equations  of  the  device  are  cast  in  the  following  form 

i 

0 


fi(v,y,v,y,u) 

h(vj,v,y,u) 


where  v  is  the  vector  of  terminal  voltages,  i  the  vector  of 
terminal  currents,  y  the  vector  of  internal  state  variables  and 
u  the  vector  of  independent  control  variables. 

A  linear  form  of  Eq.  (2)  is  used  by  VTB,  as  given  below 


i(t) 

Gn(t)v(t)  +  Gn(t)y(t)  —  B\(t) 

0 

_G2\(t)v(t)  +  G22{t)y(t)  —  B2(t)  _ 

(3) 


where  G  is  the  conductance  matrix  and  B  the  current  source 
vector.  For  a  nonlinear  device,  Eq.  (3)  is  obtained  by 
linearizing  the  nonlinear  governing  equations  with  the  Tay¬ 
lor  series  expansion. 

The  RC  modeling  approach  essentially  treats  each  electric 
device  as  a  system  of  conductors  and  current  sources.  A  RC 
model  accepts  the  values  of  terminal  voltages  as  input  and 
calculates  the  conductance  matrix  and  the  current  source 
vector  as  output.  A  dynamic  RC  model  also  needs  to  conduct 
the  integration  of  the  differential  equations  for  a  specified 
time  step.  Each  RC  model  only  handles  its  own  governing 
equations.  The  RC  numerical  solver  of  VTB  will  take  care  of 
the  solving  of  all  interconnected  devices  in  a  system.  For  a 
dynamic  system,  at  each  time  step  the  RC  solver  adjusts  the 
voltages  of  all  connecting  nodes  in  the  electric  circuit  based 
on  the  previous  values  until  the  numerical  convergence  is 
achieved,  i.e.  the  sum  of  incoming  currents  is  approximately 
zero  at  all  connecting  nodes.  A  summary  of  the  algorithm 
used  by  the  RC  solver  is  provided  in  the  Appendix  A,  which 
is  essentially  a  nonlinear  equation  solving  procedure  based 
on  Newton’s  algorithm.  More  details  can  be  found  in  the 
studies  by  Cokkinides  and  Beker  [3]. 


i  O  vi 


0  O  v0 


V] 


Fig.  4.  Diagram  of  a  battery  RC  model  where  vceii  =  v’i  —  v’o  and 
Ceil  =  i\  =  —io-  The  matrices  b  and  g  are  given  in  Eqs.  (5)  and  (6). 


source  device  with  two  terminals,  as  shown  in  Fig.  4.  For  a 
two-terminal  electric  device,  RC  modeling  expects  the 
following  relation 

i  —  gv  —  b  (4) 


The  conductance  is  given  by 

d  i 
dv 

The  current  source  is  calculated  by 


(5) 


b  —  gv  —  i 


(6) 


Put  the  relations  among  dependent  variables  into  RC 
matrix  form,  we  have 


/  =  GV  —  B 


(7) 


where 


V  = 

1  = 

G  — 

B  = 


vo 

_vi 

k 

h 


g  ~g 
~g  g 

—b 

b 


(8) 

(9) 

(10) 

(11) 


The  behavior  of  a  battery  is  usually  described  in  a  set  of 
nonlinear  equations  (for  dynamic  equations,  discretization 
in  time  domain  is  needed)  as  given  below 

/(v,  i,y)  =  0  (12) 

where y  is  a  vector  of  state  variables.  To  provide  a  RC  battery 
model,  the  conductance  matrix  and  the  current  source  vector 
must  be  obtained  from  the  governing  equations  of  the 
battery.  Two  battery  models  are  presented  subsequently  to 
demonstrate  the  details  of  the  RC  modeling  approach. 


2.  RC  modeling  of  batteries 


2.1.  A  simple  RC  model  for  general  batteries 


Battery  is  a  typical  electric  device  with  the  input/output  of 
current  and  voltage.  In  fact,  it  can  be  treated  as  a  nonlinear 


A  simple  battery  model  was  developed  as  an  extension  of 
the  widely  used  model  Eq.  (1).  The  governing  equations  of 
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the  model  is  given  by 

V  =  Eeq  +  *Rint  (13) 

where  the  equilibrium  potential  of  the  battery  is  a  function  of 
the  SOC  of  the  battery  and  can  be  fitted  directly  from  the 
experimental  data 

Eeq  =/(  SOC)  (14) 

If  the  following  assumptions  are  valid:  (i)  the  Nernst 
equation  is  a  valid  description  of  the  equilibrium  potential 
of  the  relevant  electrochemical  reactions;  (ii)  the  main 
electrochemical  reactions  on  positive  and  negative  electro¬ 
des  have  fast  kinetics;  (iii)  the  capacity  of  the  positive 
electrode  roughly  matches  that  of  the  negative  electrode 
(i.e.  no  precharge  and  surplus  capacity  on  either  electrode), 
we  have 


2.2.  RC  model  of  a  nickel-metal  hydride  cell 

The  model  presented  here  uses  the  planar  electrode 
approximation  for  both  positive  and  negative  electrodes, 
i.e.  each  electrode  is  treated  as  a  huge  planar  electrode  with  a 
surface  area  equal  to  the  total  internal  surface  area  of  the 
porous  electrode.  The  diffusion  process  inside  the  solid 
active  material  has  been  neglected.  Since  battery  electrodes 
are  normally  made  to  be  very  thin  to  minimize  ohmic 
resistance,  a  planar  electrode  model  usually  gives  insignif¬ 
icant  errors  and  costs  much  less  computation  compared  with 
a  rigorous  porous  electrode  model. 

Four  reactions  in  a  nickel-metal  hydride  cell  have  been 
considered  (the  hydrogen  oxidation  reaction  on  the  nickel 
electrode  has  been  neglected),  as  listed  below 

•  Nickel  electrode 


RT 

Eeq  —  Eq  ~\  —  In 

nF 


SOC 


1  -  SOC 


(15) 


O  NiOOH  +  H20+e“ 


discharge 

charge 


Ni(OH)2  +  OH" 


(22) 


Assuming  no  side  reactions,  then  the  relationship  between 
SOC  and  charge/discharge  rate  is  given  by 


dSOC  _  i 
d t  ~  3600Cbat 


(16) 


The  equivalent  conductance  and  current  source  of  the 
model  are  calculated  as  follows 


o  20H-  ->  i02+H20  +  2e~ 

•  Metal  hydride  electrode 

discharge 

o  MH  +  OH"  H20  +  M  +  e" 

charge 

o  j02  +  H20  +  2e~  — »  20H~ 


(23) 


(24) 

(25) 


di  _  1 
dv  Rint 


b  =  gv  —  i  = 


(17) 

(18) 


In  matrix  form,  we  have  the  RC  model  for  general 
batteries. 


/  =  GV  —  B 


(19) 


The  main  reaction  on  the  positive  electrode  is  the  redox 
reaction  of  nickel  active  material  Eq.  (22)  while  that  on  the 
negative  electrode  is  the  redox  reaction  of  metal  hydride 
material  Eq.  (24).  The  side  reaction  on  the  positive  electrode 
is  oxygen  evolution  Eq.  (23)  while  that  on  the  metal  hydride 
electrode  is  oxygen  reduction  Eq.  (25). 

The  kinetics  of  reactions  Eqs.  (22)-(24)  can  be  described 
with  the  Butler- Volmer  equation 


where 


G  — 


B  = 


-  1  1  - 

Rint  Rint 

1  1 

Rint  Rint  - 

~  Eeq~ 

Rint 
Re q 

-  Rint  - 


(20) 


(21) 


The  above  battery  model  is  actually  a  linear  model 
between  current  and  voltage.  Due  to  its  simplicity,  computa¬ 
tional  requirement  of  the  model  is  trivial.  This  model  is 
more  realistic  than  Eq.  (1)  and  can  give  reasonable  predic¬ 
tions  for  many  battery  systems.  However,  it  is  still  over¬ 
simplified.  First-principle  modeling  is  needed  to  obtain 
more  accurate  predictions  and  get  mechanistic  information 
of  a  battery.  An  example  is  presented  below  for  a  nickel- 
metal  hydride  cell. 


7  i  —  *o,i 


72  —  *0,2 


(A^pos  Req,2,ref) 


F 

RT 


(A^pos  Req,2,ref)^ 


(27) 


73  —  *0,3 


CMH 

T'MH,  ref , 


0.5E  /  a  i  0 

eXPl  -^r(A0neg-  ^,3, ref) 


exp 


Te,ref  J  \  RT 

0.5F 


RT 


(^^*neg  Req,3,ref) 


(28) 
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where  the  potential  difference  at  the  solid-liquid  interface 
on  each  electrode  is 


A0pOs  ^pos,s  ^pos,l  (29) 

A^neg  ^neg,s  ^neg,l  (20) 


Due  to  the  huge  electrochemical  driving  force  for  oxygen 
reduction  on  the  negative  electrode,  a  limiting-current  equa¬ 
tion  is  used  for  the  rate  of  reaction  Eq.  (25) 


.  _  P  o2  . 

J  4  —  M-,  ref 

P  02,ref 


(31) 


The  charge  balances  on  the  electrodes  is  described  by 
Cell  ^pos^pos^pos (y'l  T  jl)  (32) 

Cell  GiegCegAneg  (.73  T 74  )  (33) 

The  mass  balances  of  nickel  active  material  is  given  by 


d  c 


Ni(OH)2 


n 

F 


pos  d  t 

where  the  effective  thickness  of  nickel  active  material  is 


(34) 


ly,  pos  — 


^Ni(OH)2 


(35) 


PNi(OH)2^pos^pos 

The  mass  balance  of  metal  hydride  material  is  given  by 


CMH  _  0  ^Ni(OH)2^pos4^MH  ^  ^ 

^MH.max  ^yMHy4neg4TNi(OH)9 


(43) 


In  Eq.  (43)  metal  hydride  material  has  been  set  to  be  five 
percent  precharged 

For  an  applied  voltage  vceii  and  an  integration  time  step 
At,  six  dependent  variables:  Ceil,  cNi(oh)2,  cM h,  Po2,  A</>pos 
and  A0neg  need  to  be  solved  with  the  model.  Applying 
the  implicit  Euler  method  on  the  model  equations 
yields 


C,pos 


C,  neg 


(new)  (°ld) 

CNi(OH)2  “  CNi(OH)2 


At 

(new)  _  (old) 

CMH  CMH 


At 


.(new) 

A 

F 


.(new) 

A _ 

F 


(44) 

(45) 


VPST)  -  P{of  _  (apos/posAPos7fw)  +aneg/„egAneg7fw)) 


RT 

At 

F 

(46) 

.(new)  _ 

Cell  — 

a  1  A  (  Cnew)  4-  CnewC 
Ciposlpos-ripos  \  J 2  J 

(47) 

.(new) 

Cell 

—  n  1  A  (  Cnew)  4-  7-(new)\ 

Uneg'neg^neg  '  J 4  J 

(48) 

(new) 

Cell  = 

-  A<„esw)  -  +  47* 

R\ni 

(49) 

Vneg 


dCMH 
d  t 


h 

F 


(36) 


where  the  effective  thickness  of  metal  hydride  active  mate¬ 
rial  is 


C,neg  — 


Tmh 


Pmh^ 


neg^neg 


(37) 


The  mass  balance  equation  of  oxygen  is  given  by 


^gas  4/?o2 
RT  d t 


{OpOSlpOSApOS  j 2  T  ^neg^negAnegJd) 


F 


(38) 


The  separator  is  treated  as  an  ohmic  resistance,  which  gives 
Cell^int  0pos,l  ^neg,l  (39) 

The  cell  potential  is  calculated  by 


Cell  0pos,s  */Ceg,s  A0pos  A(f)ne„  T  Cell^int  (40) 


Due  to  the  comparatively  excess  amount  of  metal  hydride 
material,  the  SOC  of  a  nickel-metal  hydride  cell  is  given  by 
the  charged  state  of  nickel  active  material 


soc=i-  CNl(0H)2 

CMi(OH)2,max 


(41) 


The  initial  setting  of  SOC  for  the  model  will  be  converted 
into  corresponding  concentration  of  nickel  active  material 
and  metal  hydride  material,  i.e. 


N‘(QH)2  =1-SOC0 

cNi(OH)2  ,max 


(42) 


where  the  superscript  ‘(old)’  indicates  the  previous  values 
and  superscript  ‘(new)’  the  current  values.  It  is  worth  noting 
that  the  old  values  of  algebraic  variables  fceii ,  A</>pos  and 
A0neg  have  no  effect  on  the  solution  of  the  current  values  of 
dependent  variables  because  no  time  derivatives  of  these 
variables  appear  in  the  governing  equations.  The  above 
equations,  which  can  be  written  in  the  form  of  Eq.  (12), 
are  solved  with  the  nonlinear  equation  solver,  general  non¬ 
linear  equation  solver  (GNES)  [4]. 

Analytical  solutions  for  the  RC  matrices  are  not  available 
for  the  above  model.  The  equivalent  conductance  g  and 
current  source  b  of  the  RC  battery  model  are  obtained  by 
using  finite  difference  approximations.  That  is,  g  is  calcu¬ 
lated  by  a  finite  difference  approximation  for  di/dv  in 
Eq.  (5).  This  approximation  is  obtained  by  first  setting 
vCeii  followed  by  calculating  /ceii  and  the  other  dependent 
variables  listed  above  by  solving  Eqs.  (44)-(49).  Next,  a 
small  voltage  Avcen  is  added  to  vceu  and  the  calculation  is 
repeated  to  obtain  Ceil  +  A  Ceil-  These  values  are  used  to  give 
g  and  b 


8  = 


(Cell  T  A  Cell )  Cell 

(  Vceii  T"  A Vcen)  Vceu 


b  —  g^cell  Cell  — 


(Cell  T  A/Cell)  Cell 
(Vceii  T  A Vceii)  VCell 


Vceii  Cell 


(50) 

(51) 


In  matrix  form,  we  have  the  RC  model  of  the  nickel-metal 
hydride  cell  as  follows 


/  =  GV  —  B 


(52) 
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where 


G  — 


B  = 


(kell  +  Ai'cell)  kell  (^'cell  +  A/Cep)  ^'cell 

(l^cell  +  AvCell)  VCell  ( ^ cel  1  +  Avcep)  ^cell 

O'cell  T  A/cep)  kell  (kell  T  A/cep)  kell 

(+ell  +  Avcell)  Vcen  (vCeii  +  Avcep)  Vcen  _ 

(53) 


(kell  +  Afceii)  ^cell 

7  "  7  7  l^cell  +  icell 

(+cell  +  AVceiiJ  Vcell 

(kell  +  A/ceii)  /ceu 

+ell  kell 


l^cell  3“  Avcen)  Vceu 


(54) 


The  perturbation  of  the  cell  voltage  Avceii  must  be  small 
enough  (e.g.  1  x  10”4  V)  to  provide  accurate  numerical 
approximations  for  G  and  B.  The  numerical  approximation 
of  RC  matrices  essentially  allows  more  complex  battery 
models  to  be  treated  as  a  linear  model  between  current  and 
voltage  at  a  specific  time  point. 

Good  efficiency  and  robustness  are  highly  desirable  for  a 
first-principle  RC  battery  model.  It  is  obvious  that  even  with 
the  planar  electrode  approximation,  the  first-principle  model 
of  a  nickel-metal  hydride  cell  is  highly  nonlinear  and 
nontrivial  to  solve.  In  the  simulation  of  a  whole  power 
system,  since  a  battery  model  may  be  called  hundreds  or 
even  thousand  times  by  the  RC  solver,  a  slow  model  could 
cause  a  simulation  to  last  several  days,  which  will  be 
intolerable  for  most  projects.  Normally,  simplified  models 
and  good  numerical  solvers  should  be  used  to  achieve 
satisfactory  calculation  efficiency.  In  the  RC  approach, 
the  input  to  a  battery  models  is  terminal  voltage  values; 


however,  the  behavior  of  a  battery  is  highly  sensitive  to  the 
applied  voltages,  as  is  obvious  from  the  exponential  depen¬ 
dence  of  battery  current  on  the  voltage  values  in  Eqs.  (26)- 
(28).  Some  guess  values  from  the  RC  solver  can  easily  lead 
to  numerical  difficulty  of  a  battery  model.  Thus  a  RC  battery 
model  needs  to  be  built  robust  enough  to  avoid  numerical 
failures  of  the  whole  simulation. 


3.  Results  and  discussions 

To  incorporate  a  RC  model  into  VTB,  a  dynamic  link 
library  (DLL)  hie  needs  to  be  constructed  based  on  the  C++ 
interface  specified  by  Cokkinides  and  Beker  [3].  The  DLL 
hie  of  a  RC  model  can  be  debugged  independently.  VTB  will 
accept  and  use  the  DLL  hie  of  a  user-built  RC  model  as  it 
does  with  native  RC  models.  Communication  between  a  RC 
model  and  the  VTB  is  through  the  exchange  of  model 
parameters  and  RC  matrices.  Sample  source  code  has  been 
provided  by  Cokkinides  and  Beker  [3]  to  facilitate  the 
construction  of  RC  models  for  the  VTB  environment.  A 
C++  class  of  a  battery  model  is  derived  from  the  base  class 
CCompanionDevice.  The  following  subroutines  need  to  be 
dehned  (supposing  BatModel  is  the  name  of  our  battery 
model  class): 

•  BatModel:  :BatModel(),  defines  model  parameters  and 
states; 

•  BatModel:  :GuiToInternal(),  converts  user  settings  to 
internal  model  parameters; 

•  BatModel:: Get VectorIcon(short*  vicondata),  provides  the 
icon  of  the  model; 


External  Device  Properties 


El 


Name 
r  Parameters 


DLL 


Battery 


□ 


Capacity 

1.0000 

A-hr 

Soc 

0.25000 

Nominal  Potential 

1.2500 

V 

Electron  T ransfer  Number 

1.0000 

Internal  Resistance 

10.000 

milliOhm 

J 


1 

1  3 

Nodes 

Cancel 

Comment 

General  battery  model 

OK 

Fig.  5.  The  parameter  setting  dialog  for  the  simple  RC  model  of  general  batteries. 
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Table  1 

Parameters  for  the  simplified  model  of  general  batteries 

Parameter 

Value 

C>at  (A  h) 

1.0 

£o(V) 

1.25 

n 

1.0 

Rmt  (Q) 

1.0  x  10~2 

Battery  1  Battery2 


Fig.  6.  Worksheet  of  sample  simulation  I:  two  batteries  with  initial  SOC  of 
25%  connected  with  a  battery  with  initial  SOC  of  85%. 


•  B atModel : : TdsInitExt(DE V_INTERFACE*  tdsi,  double* 
vol,  double*  beq,  double**  con),  initializes  model  states 
and  provides  VTB  with  initial  RC  matrices; 

•  BatModel::TdsStepExt(DEV_INTERFACE*  tdsi),  calcu¬ 
lates  the  model  states  and  RC  matrices  based  on  the 
voltage  input  from  VTB  RC  solver. 

After  the  DLL  of  a  RC  model  has  been  developed,  it  needs 
to  be  put  into  a  specified  directory  to  be  accepted  by  VTB. 
Using  a  RC  model  is  straightforward  in  VTB.  To  create  a 
model  object,  one  selects  the  corresponding  model  item 
from  menus  or  toolbars.  The  model  object  is  represented  by 
an  icon  and  can  be  moved  around  in  the  VTB  worksheet.  The 


Batteryl 


Fig.  8.  Worksheet  of  sample  simulation  II:  a  capacitor  and  an  inductor 
connected  in  series  with  a  battery  with  initial  SOC  of  95%. 


parameters  of  an  object  can  be  easily  modified  through  a 
dialog  box  (see  Fig.  5)  invoked  by  double-clicking  the  icon 
of  the  object. 

Two  simulations  are  given  for  the  simple  RC  model  of 
general  batteries.  The  parameters  of  the  battery  model  are 
listed  in  Table  1.  In  Fig.  6,  three  batteries  are  connected 
together.  The  SOCs  of  the  first  two  batteries  are  set  to  be 
25%,  and  the  SOC  of  the  third  one  is  set  to  be  85%.  However, 
as  shown  in  Fig.  7,  it  is  predicted  that  two  less  charged 
batteries  connected  in  series  can  actually  produce  a  voltage 
high  enough  to  charge  a  near  fully  charged  battery  in  15  s, 
which  is  quite  an  unintuitive  result  (and  a  fact  in  reality).  In 
Fig.  8,  a  battery  with  an  initial  SOC  of  95%  is  connected 
with  a  1C  capacitor  and  a  1H  inductor.  The  simulation 
results  in  Figs.  9  and  10  show  fast  transient  behavior  of 
voltage  and  current.  The  capacitor  is  charged  while  the 
battery  is  discharged  at  first.  When  the  capacitor  is  near 
fully  charged,  the  current  in  the  circuit  approaches  zero  (the 
capacitor  is  essentially  open-circuit  after  fully  charged)  and 
the  battery  voltage  approaches  the  open  circuit  potential  of 
the  battery.  Since  the  battery  discharge  process  only  lasts  for 
10  s,  the  SOC  of  the  battery  only  slightly  decreases,  as 
shown  in  Fig.  11. 


Fig.  7.  Predicted  voltage  of  battery  3  in  simulation  I. 
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Table  2 


Parameters  for  the  model  of  a  nickel-metal  hydride  cell 


Parameter 

Value 

apos  (cm2/cm3) 

4000.0 

aneg  (cm2/cm3) 

1000.0 

Ap0s  (cm-) 

325.0 

Aieg  (cm-) 

360.0 

ce  (mol/cm3) 

7.0 

X 

10~3 

ce,ref  (mol/cm3) 

1.0 

X 

10“3 

CMH,max  (mol/cnT) 

1.0 

X 

10”1 

CMH,ref  (CMH.max) 

0.5 

CNi(OH),,max  (mol/cm  ) 

3.7 

X 

10~2 

cNi(OH)2,ref  (cNi(OH)2,max) 

0.5 

/„  |  (A/cm2) 

1.0 

X 

10~4 

i0,2  (A/cm2) 

1.0 

X 

10”10 

*0,3  (A/cm2) 

bo 

X 

10~4 

k.,4  (A/cm2) 

1.0 

X 

10”4 

/pos  (cm) 

3.3 

X 

10~2 

/neg  (cm) 

,  .  Os 

2.8 

X 

10”2 

- O 

^ni(oh)2  (g/cm2)  6.8  x  10  2 

LMh  (g/cm2)  1.13  x  10  1 


A/Ni(OH)2  (g/mol) 

92.71 

Mmh  (g/mol) 

70.58 

P°o2  (atm) 

1.0  x  10"2 

PO2  ,ref  (atm) 

1.0 

Rint  («) 

1.0  x  10"3 

SOC0 

0.95 

At  (s) 

1.0  x  10_1 

T( K) 

298.15 

t/eq,l,ref  (V) 

0.527 

t/eq,2,ref  (V) 

0.4011 

t/eq,3,ref  (V) 

-0.8279 

Vgas  (cm3) 

PNi(OH)2  (g/cm3) 

Pmh  (g/cm3) 

1.0  x  10_1 
3.4 

7.49 

NiMHI 


100.00m 


Fig.  13.  Worksheet  of  sample  simulation  III:  a  nickel-metal  hydride  cell 
with  initial  SOC  of  95%  discharges  through  a  100  m  Q  resistor. 


Two  simulations  are  provided  for  the  RC  model  of  a 
nickel-metal  hydride  cell.  The  parameters  of  the  nickel- 
metal  hydride  cell  model  are  listed  in  Table  2.  The  para¬ 
meters  used  as  the  ‘old’  variables  for  the  initial  conditions 
are  shown  with  a  superscript  ‘O’  in  Table  2.  Fig.  12  shows  the 
parameter  dialog  for  the  model  in  VTB.  In  Fig.  13,  a  nickel- 
metal  hydride  cell  with  an  initial  SOC  of  95%  is  connected 
with  a  constant  resistor  of  100  m  Q.  The  simulated  current  in 
circuit  is  shown  in  Fig.  14.  It  is  observed  that  with  the  drain 
of  SOC  of  the  nickel-metal  hydride  cell,  the  current  in  the 
circuit  gradually  decreases,  and  at  the  end  of  the  discharge 
process,  the  current  drops  rapidly  to  zero.  Since  the  cell  is 
connected  to  a  constant  resistor,  the  cell  voltage  curve 
similar  to  Fig.  14  can  be  expected. 

The  simulated  SOC  and  hydrogen  concentration  in  the 
metal  hydride  electrode  are  shown  in  Figs.  15  and  16, 
respectively.  Approximately  constant  depletion  rate  of 
the  cell  capacity  is  predicted.  It  is  shown  that  only  half 
of  metal  hydride  active  material  is  utilized,  which  is 


Parameters 


Initial  SOC  (State  of  Charge) 

0.95000 

▲ 

Ambient  temperature 

25.000 

Degree  C 

Rated  capacity 

4.5000 

A.h 

Internal  electrical  resistance 

1.0000 

milliOhm 

K0H  molar  concentration 

7.0000 

M 

Geometry  area  of  the  nickel  electrode 

325.00 

cmA2 

Thickness  of  the  nickel  electrode 

33.000 

millicm 

Loading  of  nickel  active  material 

G8.000 

millig/crrT2 

Specific  surface  area  (Pos) 

4.0000 

kilocmA2/crrF3 

Geometry  area  of  the  metal  hydride  electrode 

360.00 

cm  2 

Thickness  of  the  metal  hydride  electrode 

28.000 

millicm 

Loading  of  metal  hydride  active  material 

0.11300 

g/crrF2 

d 

Specific  surface  area  (Nea) 

1.0000 

kilpcmA2/cmA3 

1  3 

Nodes  I  -i  .-| 

Cancel 

Comment  Nickel/metal-hydride  battery  model 

OK 

Fig.  12.  The  parameter  setting  dialog  for  the  nickel-metal  hydride  cell  model. 


Time  (s) 

Fig.  14.  Predicted  current  of  the  cell  in  simulation  III. 


Time  (s) 

Fig.  15.  Predicted  SOC  of  the  cell  in  simulation  III. 


Time  (s) 


Fig.  16.  Predicted  mol  fraction  of  hydrogen  on  the  metal  hydride  electrode  of  the  cell  in  simulation  III. 
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0  500  1000  1500  2000  2500  3000 

Time  (s) 

Fig.  17.  Predicted  oxygen  pressure  of  the  cell  in  simulation  III. 


consistent  with  the  positive-limited  design  of  the  battery. 
The  predicted  oxygen  pressure  in  the  nickel-metal 
hydride  cell  is  given  in  Fig.  17.  It  indicates  that  the  oxy¬ 
gen  pressure  decreases  gradually  with  discharge  time. 
As  shown  in  Fig.  18,  the  potential  difference  at  the  solid- 
liquid  interface  on  the  positive  nickel  electrode  gradually 
decreases  at  first,  and  at  the  end  of  the  discharge  process, 
it  changes  rapidly  to  a  negative  value.  In  Fig.  19,  however, 
it  is  shown  that  the  potential  difference  at  the  solid- 
liquid  interface  on  the  negative  metal  hydride  electrode 
only  slightly  changes  during  the  whole  discharge  process. 
It  is  observed  that  at  the  end  of  the  discharge  process, 
the  potential  difference  at  the  solid-liquid  interface  on 


the  positive  electrode  is  equal  to  that  on  the  negative 
electrode  and  causes  zero  current  across  the  cell  (and  whole 
circuit). 

A  simulation  involving  two  nickel-metal  hydride  cells  is 
shown  in  Fig.  20.  The  initial  SOCs  of  two  cells  are  set  to  be 
different.  The  initial  SOC  of  the  first  one  is  95%  and  that  of 
the  second  one  is  75%.  The  predicted  discharge  current  is 
shown  in  Fig.  21.  As  can  be  expected,  due  to  the  high  voltage 
of  two  serially  connected  nickel-metal  hydride  cells,  the 
current  in  the  circuit  almost  doubles  that  with  a  single 
nickel-metal  hydride  cell.  The  predicted  SOCs  in  two  cells 
are  given  in  Figs.  22  and  23.  It  is  interesting  to  note  that  the 
less  charged  cell  determines  the  utilization  of  both  cells,  i.e. 


Fig.  18.  Predicted  potential  difference  at  solid-liquid  interface  on  the  nickel  electrode  of  the  cell  in  simulation  III. 
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Time  (s) 

Fig.  19.  Predicted  potential  difference  at  solid-liquid  interface  on  the  metal  hydride  electrode  of  the  cell  in  simulation  III. 


NiMHI  NiMH2 


Fig.  20.  Worksheet  of  sample  simulation  IV:  two  nickel-metal  hydride 
cells  with  different  initial  SOCs  discharge  through  a  100  m  Q  resistor. 

when  the  second  nickel-metal  hydride  cell  gets  fully  dis¬ 
charged,  it  blocks  the  discharge  of  the  unused  capacity  of  the 
first  cell.  Thus,  it  confirms  a  common  sense  in  the  battery 


industry,  i.e.  it  is  not  wise  to  use  two  cells  with  different 
SOCs. 

RC  modeling  easily  allows  simulations  with  multiple 
objects  of  the  same  model,  i.e.  once  a  RC  model  is  built, 
it  can  be  repeatedly  used  for  many  objects  in  a  complex 
system  model.  This  is  somewhat  different  from  the  conven¬ 
tional  battery  modeling  approach,  e.g.  if  two  connected 
batteries  are  to  be  simulated,  a  more  complex  model  based 
on  the  single  battery  model  is  required,  which  is  usually  a 
nontrivial  task.  Usually  different  objects  of  the  same  model 
may  have  different  parameters,  and  they  are  handled  simul¬ 
taneously  with  the  RC  solver.  It  is  obvious  that  challenging 
battery  simulations  can  be  more  easily  conducted  with 
the  RC  modeling  approach  in  a  simulation  environment 


Fig.  21.  Predicted  current  of  the  cell  in  simulation  IV. 


B.  Wu  et  al.  /  Journal  of  Power  Sources  93  (2001)  186-200 


199 


Fig.  22.  Predicted  SOC  of  the  first  nickel-metal  hydride  cell  in  simulation  IV. 


Fig.  23.  Predicted  SOC  of  the  second  nickel-metal  hydride  cell  in  simulation  IV. 


Fig.  24.  Schematic  diagram  of  the  RC  model  objects  for  the  worksheet  in 
Fig.  20. 


(e.g.  in  VTB)  than  with  the  direct  modeling  and  coding 
approach. 


4.  Conclusions 

The  simulation  of  interactions  among  batteries  and  other 
electric  devices  is  usually  a  nontrivial  task.  In  this  study,  the 
RC  battery  modeling  for  electric  circuit  simulation  is  pre¬ 
sented.  Two  RC  battery  models  (a  general  battery  model  and 
a  nickel-metal  hydride  cell  model)  are  demonstrated  on  the 
VTB  platform.  It  is  shown  that  challenging  simulations 
involving  batteries  and  other  electric  devices  can  be  easily 
handled  with  the  RC  modeling  approach. 
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Appendix  A.  Algorithm  of  the  RC  solver  in  VTB 


For  an  electric  circuit  system,  from  the  charge  conserva¬ 
tion  and  RC  model  equations,  we  have  the  governing 
equations  for  the  dependent  voltage  variables  v  at  connecting 
nodes. 

n 

F(v)  =  +b)  =  0  (A-1) 

i=i 


The  Jacobian  matrix  of  the  above  equation  with  respect  to 
v  is  given  by 


J(0 


(A.2) 


The  Newton  iteration  can  be  used  to  solve  the  dependent 
voltage  variables. 

vk+i  =vk  _j\  (v*)  (A.3) 

To  get  the  new  iteration  values  with  Eq.  (A.3),  the 
inversion  of  Jacobian  matrix  is  not  conducted  due  to  the 
efficiency  and  robustness  considerations.  Instead,  LU 
decomposition  procedure  is  used.  First,  J(pck)  is  decomposed 
into  upper  and  lower  triangle  matrices 

J{vk)  =  LU  (A.4) 

Then  the  following  equation  is  obtained 

LU(vk  -  vM)  =F(vk)  (A.5) 

which  can  be  solved  by  forward  and  backward  substitution 
processes,  i.e.  solve 

Lm  =  F(vk)  (A.6) 


for  vector  m,  then  solve 


U(vk-vk+1)=m  (A.7) 

for  vk  —  vi+ 1  and  finally  get  vi  + 1 . 

For  example,  to  solve  the  circuit  of  three  RC  model 
objects  connected  together  in  Fig.  20,  we  have  three  depen¬ 
dent  variables  vo,  Vi  and  V2  (see  Fig.  24).  The  governing 
equations  for  them  are 

40)  +  4°  =  0  (A.8) 
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Substitution  of  Eqs.  (A.ll)-(A.13)intoEqs.  (A.8)-(A.10) 
yields 
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It  is  obvious  that  as  long  as  the  RC  parameters  are 
provided  by  RC  models,  the  values  of  vo,  vi  and  V2  can 
be  easily  solved  from  Eq.  (A.  14).  More  complex  electric 
circuits  can  be  solved  in  the  same  way.  The  construction  of 
the  governing  equations  like  Eq.  (A.  14)  of  the  whole  simu¬ 
lated  system  is  automatically  handled  by  VTB. 
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